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基于 多 源 遥 感 数据 的 1970 一 2020 年 巴尔 喀什 湖 动 态 监测 
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摘 要 : 利用 ICESat-1、CryoSat-2 测 高 数据 与 水 文 站 实测 数据 ,分 析 了 1970 一 2020 年 巴尔 喀什 湖水 


位 变化 特征 , 结 


合 MOD09A1 数据 及 相关 资料 分 析 了 巴尔 喀什 湖 长 时 间 序 列 酝 
并 利用 巴尔 喀 n 湖 流域 灌溉 数据 及 Climatic Research Unit(CRU) 气 象 数 据 ,简要 探究 


只 和 水 量 变化 状况 ， 
究 了 气候 变化 及 


人 类 活动 对 湖泊 长 期 变化 的 影响 。 结 果 表 明 : 巴 尔 喀 什 湖 水 位 在 研究 时 段 内 总 体 呈 现 先 下 降 后 波 


动 上 升 的 趋势 ,1987 年 水 位 值 达 到 最 低 (340.64 m), 年 内 水 位 增长 主要 发 生 在 2 月 底 至 
季 (4 一 10 月) 水 位 波动 比 冷 季 (11 月 一 次 年 3 月 ) 更 剧烈 ;1970 一 2020 年 巴尔 喀什 湖 郑 
缩小 约 16.77%; 在 此 期 间 巴 尔 喀 什 湖水 量 损失 约 为 12.33 kmo IA 


km? J 2? F| 16641.93 km’, X if 
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萨克斯 坦 境内 灌区 耕地 扩张 等 人 类 活动 引起 ，; 气温 与 降水 在 总 体 趋势 上 与 湖泊 的 动态 变化 并 无 明 
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湖泊 对 气候 与 环境 变化 极为 敏感 ,在 维持 区 
域 生 态 平 衡 中 起 着 重要 作用 ”, 湖 泊 的 水 位 、 面 积 和 
水 量变 化 与 区 域 气温 、 降 水 .蒸发 及 人 类 活动 等 因 
素 密切 相关 ,因而 其 动态 变化 可 间接 反映 区 域 气 
候 变化 以 及 人 类 活动 对 湖泊 生态 环境 的 影响 “3 
巴尔 喀什 湖 是 中 亚 干 旱 区 第 三 大 水 体 ， 20 世 纪 80、 
90 年 代 , 巴 尔 喀什 湖 流域 内 不 合理 的 人 类 活动 使 其 
生态 环境 遭 到 严重 破坏 ,湖水 水 位 大 幅 下 降 湖泊 
面积 萎缩 、 土 地 生产 力 下 降 “。2008 年 ,联合 国 开发 
计划 署 在 地 理 年 鉴 中 表明 巴尔 喀什 湖面 临 着 完全 
于 酒 的 可 能 性 ”, 随 后 该 湖泊 的 生态 系统 持续 吸引 
着 公众 的 关注 ,因而 对 巴尔 喀什 湖 及 其 环境 的 长 期 
监测 至 关 重 要 。 

湖泊 动态 变化 主要 体现 在 湖泊 面积 和 水 位 的 
变化 。 传 统 的 湖泊 水 位 监测 依靠 水 文 站 点 实测 ， 
这 种 方法 虽然 精度 较 高 ,但 对 人 力 和 物力 的 消耗 
较 大 ” , 且 受 限于 自然 条 件 与 空间 分 布 等 因素 ,地 域 
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电量 讨论 各 因素 对 湖泊 动态 变化 的 影响 ,; 
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还 需 进 一 步 研 究 。 


偏远 的 内 陆 湖泊 难以 获得 长 时 间 的 连续 观测 资 
料 。 此 外 ,传统 水 位 监测 共享 程度 较 低 ,增加 了 数 
据 获取 的 难度 ”。 卫 星 测 高 技术 作为 1 种 空间 大 地 
测量 技术 , 现 被 广泛 应 用 于 内 陆 湖 泊 的 水 位 监测 及 
分 析 其 对 环境 的 响应 ,与 传统 的 水 位 监测 方式 相 
比 ,卫星 测 高 技术 能 够 大 范围 .周期 性 地 监测 陆地 
水 体 的 各 种 动态 变化 "" ,多 项 研究 结果 表明 卫星 测 
高 技术 能 够 有 效 监测 CS eee k 卫 
星 测 高 数据 主要 来 源 于 激光 与 雷达 2 种 类 型 的 高 度 
计 。ICESat-1 是 世界 上 第 一 颗 星 载 激光 测 高 卫星 ， 
它 利 用 地 球 科 学 激光 测 高 系统 测量 全 球 地 表 高 
度 , 其 地 面 足 迹 点 的 直径 约 为 70 m, 沿 轨迹 相 邻 足 
迹 点 间距 约 为 172 m, 能 够 对 湖泊 水 位 进行 高 精度 
B5 dj 59". CryoSat-2 卫星 由 ESA (European Space 
Agency) 于 2010 年 发 射 , 拱 载 的 是 合成 孔 色 干涉 雷 
达 高 度 计 SIRAL, 沿 轨迹 相 邻 足迹 点 间距 约 为 300 
m, 赤 道上 的 相 邻 轨迹 间距 约 为 8 km ,与 传统 的 雷达 
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测 高 系统 如 ERS-L2 Envisat GFO VA & Jason- 1/2 £F 
相 比 ,CryoSat-2 轨道 更 密集 ,是 迹 点 间距 与 沿 轨 间 
距 更 小 ,可 以 监测 更 多 的 小 型 内 陆 湖 泊 "”。ICESat- 
1 的 运行 时 间 为 2003 一 2009 4E. , CryoSat-2 的 运行 时 
间 为 2010 年 至 今 , 因 而 有 多 项 研究 结合 这 两 颗 卫 星 
的 观测 数据 构建 了 长 时 间 序 列 的 观测 资料 。 宋 春 
桥 等 "中 利用 ICESat-1 与 CryoSat-2 卫 星 构建 了 青 
藏 高 原 131 个 湖泊 2003—2014 年 的 长 时 间 水 位 序 
列 , 且 用 ICESat-1 与 CryoSat-2 卫星 对 纳木错 湖 
2003—2014 年 的 水 位 趋势 变化 进行 了 进一步 的 分 
析 ; 黎 鹏 等 所 监测 了 2003 一 2017 年 洞庭 湖 的 水 位 变 
化 情况 。 

巴尔 喀什 湖水 位 数据 实测 记录 始 于 1938 年 ,更 
早 之 前 的 延长 资料 由 前 苏联 学 者 通过 野外 实地 考 
察 ` 水量 平衡 估算 和 相关 资料 考证 等 多 种 方式 获 
得 '” ,目前 有 很 多 研究 采用 了 巴尔 喀什 湖水 位 数 
据 。 杨 川 德 2 利 用 延长 水 位 资料 与 实测 资料 讨论 
了 1879 一 1986 年 巴尔 喀什 湖水 位 变化 及 原因 ; 龙 爱 
华 等 2 将 巴尔 喀什 湖 1879 一 2009 年 的 水 位 数据 与 
该 时 段 内 径流 量 以 及 气温 降水 等 数据 相 结合 ,对 巴 
尔 喀什 湖水 量 平衡 进行 了 人 研究 ; 郭 利 丹 等 ”讨论 了 
巴尔 喀什 湖 1810 一 2006 年 的 水 位 周期 变化 与 流域 
气候 变化 特征 的 关系 ;Propastin 778] H B. A dll 
高 资料 和 流域 气象 站 气象 资料 ,研究 了 1992—2010 
年 巴尔 喀什 湖水 位 变化 及 其 与 气候 变化 的 关系 ,并 
建立 了 湖泊 水 位 与 伊犁 河流 量 监测 的 简易 模型 。 

目前 对 于 巴尔 喀什 湖面 积 的 监测 主要 采用 遥 
FRA ANF Be. IE E E RHE TK AR EY SG C5] TE 
性 ,以 9 个 时 段 的 Landsat 影 像 为 数据 源 进行 水 体 信 
息 提取 ,讨论 了 巴尔 喀什 湖 1975 一 2014 年 面积 变化 
特征 及 其 驱动 因素 ;成 晨 等 的 采用 同 季 相 4 期 Land- 
sat 遥感 影像 ,通过 目 视 解 译 提 取 包 括 巴 尔 喀什 湖 在 
内 的 中 亚 地 区 7 个 湖泊 1978 一 2010 年 的 面积 变化 ; 
Tey EAE AE? ALF Landsat 遥感 影像 ,监测 了 巴尔 喀什 
WH 1977 一 2014 年 面积 及 湖岸 的 动态 变化 。 

以 往 对 于 巴尔 喀什 湖 的 研究 主要 以 单一 遥感 
数据 源 为 基础 ,结合 气象 数据 分 析 湖 泊 近 期 水 文 变 
化 ,对 于 面积 ,水 位 及 水 量 的 长 期 连续 变化 缺乏 较 
为 全 面 和 系统 的 分 析 , 且 目前 大 部 分 水 位 研究 都 是 
采用 水 文 站 点 的 实测 数据 或 是 基于 Jason-1、T/P 等 
高 度 计数 据 ,在 数据 精度 与 时 间 尺 度 上 都 不 够 完 
38 ,数据 源 较为 单一 ,并 未 将 传统 数据 .激光 以 及 雷 


达 测 高 数据 相 融 合 ,得 到 更 加 连续 完整 的 长 时 间 观 
测 数据 。 因 此 ,本 文 利 用 ICESat-1、CryoSat-2 卫星 
测 高 数据 及 水 文 站 实测 数据 ,构建 巴尔 喀什 湖 
1970 一 2020 年 的 水 位 序列 ,分 析 其 水 位 的 长 期 和 季 
节 变 化 特征 ,并 结合 MOD09A1 数 据 及 相关 资料 分 
析 巴 尔 喀什 湖 长 期 面积 和 水 量变 化 特征 ,同时 简要 
探究 气候 变化 及 人 类 活动 对 湖泊 动态 变化 的 驱动 
作用 ,为 巴尔 喀什 湖 的 湖泊 变化 长 期 监测 水 资源 
合理 利用 开发 及 生态 系统 保护 提供 一 定 的 科学 
依据 。 


1 研究 区 概况 


巴尔 喀什 湖 (73°21’~79°30'E, 44°45'~46°44'N ) 
位 于 哈萨克 斯 坦 东 南部 ,是 典型 的 平原 尾 间 湖 
(图 1) ,湖泊 面积 宽广 , 约 为 1.83x10' km? ,最 宽 处 约 
71 km ,总 长 约 600 km" , FAL EA ER HT 8] 
的 主要 水 动脉 经 伊犁 河 三 角 洲 注入 湖泊 西部 ,对 入 
湖 径流 量 的 贡献 率 约 为 78%2; 湖 东部 入 湖 径 流 
少 , 且 蒜 发 量 大 ,因而 东部 湖泊 矿 化 程度 明显 高 于 
Vg RE EL UU ,该 湖北 岸 水 体 很 少 , 地 势 较 高 ,多 为 
岩石 覆盖 ;南岸 小 型 湖沼 较 多 ,地 势 低 四 ,多 为 沙 
地 ”1。 伊 犁 河 一 巴尔 喀什 湖 流域 是 世界 上 最 大 的 
湖泊 生态 系统 之 一 ,流域 内 汇 人 湖泊 的 主要 河流 
有 伊犁 河 .卡拉 塔 尔 河 、 列 普 西 河 以 及 阿 亚 古 效 河 ， 
由 于 来 自 印 度 洋 ,太平洋 的 暖 湿 气流 无 法 深入 流域 
所 处 的 欧 亚 大 陆 腹地 ,因此 巴尔 喀什 湖 流域 呈现 温 
带 大 陆 性 气候 ,干燥 少雨 。 


2 数据 与 方法 


2.1 数据 来 源 

2.1.1 ICESat-1 2] à; #4 — ICESat- 1/GLAS ( Geo- 
science laser altimeter system) / 4] F 2003 4E 1 H 13 
H 7" ,通过 地 球 科学 激光 测 高 系统 (GLAS ) 回 星 下 点 
发 射 激光 信号 ,根据 激光 信号 在 卫星 与 星 下 点 之 间 
往返 所 用 的 时 间 ,来 计算 二 者 之 间 的 距离 ,进而 获 
得 星 下 点 的 高 程 ,该 卫星 获取 的 信息 时 间 范 围 为 
2003 年 2 月 20 日 至 2009 年 10 月 10 日 ,共有 640 多 
天 的 观测 数据 ™, 地 学 激光 测 高 系统 的 数据 产品 共 
有 15 种 ,本 文采 用 的 是 GLA14 全 球 测 高 数据 产品 ， 
该 数据 集中 包含 有 激光 测 高 脚 点 的 经 纬度 、 高 程 以 
及 各 项 校正 参数 等 。 
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2.1.2 CryoSat-2 测 高 数据 CryoSat-2 于 2010 年 4 
月 8 日 发 射 升 空 ,该 卫星 的 测量 模式 分 别 为 低 分 辩 
率 模式 (LRM) 、 合 成 孔径 雷达 模式 (SAR) 以 及 合成 
孔径 雷达 干涉 测量 模式 (SIRAL), CryoSat-2/SIR- 
AL 有 0 级 数据 .lb 级 数据 .2 级 数据 和 高 级 数据 4 种 
数据 产品 ,其 中 2 级 GDR 产品 包括 了 测量 时 间 、 地 
理 位 置 . 以 及 经 过 了 仪器 校正 \ 传 输 延迟 改正 、 几 何 
改正 和 地 球 物理 改正 之 后 的 高 度 信息 ,是 单独 的 全 
球 性 的 全 轨道 数据 ,即将 3 种 模式 的 测量 结果 通过 
不 同 的 处 理 过 程 ,按照 时 间 先 后 顺序 进行 整合 ,得 
到 格式 统一 的 数据 记录 。3 种 模式 的 数据 采用 了 不 
同 的 波形 重 跟 踪 算 法 来 获得 高 度 值 ,分 别 为 Re- 
fined CFI, UCL fll Refined 0COG5 ,通过 对 比 3 种 算 
法 在 巴尔 喀什 湖 湖 区 内 获得 高 度 异 常 值 的 数目 ,最 
终 采 用 CDR 产品 的 Refined OCOG 重 跟踪 算法 来 获 
取 巴 尔 喀什 湖 2010 一 2020 年 的 湖泊 水 位 变化 。 

2.1.3 湖泊 面积 数据 ”中 分 状 率 成 像 光 谱 仪 MODIS 
搭载 于 TERRA 和 AQUA 卫星 上 ,是 为 全 球 动态 测量 
提供 大 范围 数据 的 重要 传感器 ,MODIS 共 有 36 个 观 
测 波段 昌 具 有 较 高 的 分 辨 率 ,因而 被 广泛 用 于 水 体 
特征 信息 的 提取 ,目前 已 有 多 项 研究 表明 MODIS 数 
据 在 湖泊 面积 的 监测 中 具有 较 高 的 精度 。MODIS 
数据 产品 极为 丰富 ,其 中 MOD09( 来 源 于 https://lad- 
sweb.modaps.eosdis.nasa.gov/search/) 为 MODIS 陆地 2 
级 数据 产品 ,本 文选 用 的 是 MOD09 每 8 d 合 成 的 数 


注 :冰川 数据 来 源 于 中 国 第 二 次 冰川 编目 ,灌溉 数据 来 源 于 联合 国 粮 食 及 农业 组 织 。 
图 1 巴尔 喀什 湖 流域 示意 图 
Fig. 1 Location map of the Balkhash Lake Basin 


据 , 即 MOD09A1 ,空间 分 辩 率 为 5300 mo 

2.1.4 气象 数据 气象 资料 选择 了 英国 East Anglia 
大 学 Climatic Research Unit (CRU) 发 布 的 CRUTS 
v4.05 [CRU Time series version 4.05 ,获取 地 址 :CRU 
TS Version 4.05 (uea.ac.uk) |] 。 该 数据 集 是 由 CRU 收 
集 整理 多 种 数据 库 中 的 相关 资料 , 重 构 的 一 套 覆 
盖 全 球 范围 的 、 高 分 辩 率 的 连续 地 表 气 候 要 素数 
据 集 ,数据 的 空间 分 辩 率 为 0.5°x0.5°, 已 有 研究 
表明 该 数据 可 用 于 讨论 中 亚 地 区 的 气候 变化 状况 ”， 
由 于 巴尔 喀什 湖 附近 水 文 站 点 公布 的 气象 资料 不 
完整 ,因此 本 文 使 用 该 数据 集 1970 一 2019 年 的 气温 
和 降水 量 数据 作为 主要 气象 因子 来 分 析 巴 尔 喀什 
湖水 位 对 气候 变化 的 响应 (2020 年 数据 尚未 公布 )， 
同时 ,选择 了 来 源 于 NOAA 全 球 气象 站 点 数据 网 站 
(https://www.climate.gov/maps-data ) 的 部 分 站 点 数据 
来 验证 CRUTS v4.05 数据 集 在 巴尔 喀什 湖 流 域 的 适 
用 性 。 

2.4.5 其 他 数据 巴尔 喀什 湖 湖 泊 矢 量 边界 数据 来 
源 于 全 球 湖泊 和 湿地 数据 库 (Global Lakes and Wet- 
lands Database ,GLWD)1 级 数据 。 巴 尔 喀什 湖 1970 
年 起 的 实测 水 位 数据 来 源 于 HYDROLARE (http:// 
hydrolare.net/) , 其 中 1998 年 与 1999 年 数据 存在 部 
分 月 份 的 缺失 。 研 究 中 使 用 的 巴尔 喀什 湖 流域 灌 
溉 面积 来 源 于 联合 国 粮食 及 农业 组 织 FAO (http:// 


www.fao.org/aquastat/zh/countries- and- basins/regional- 
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overviews/central-asia/) 和 相关 人 研究 资料 。 
2.2 研究 方法 

本 文 的 技术 路 线 如 图 2 所 示 , 首 先 对 ICESat-1 
GLA14 和 CryoSat-2 GDR 测 高 数据 进行 预 处 理 , 提 
取 卫 星 的 足迹 点 高 程 ,筛选 出 巴尔 喀什 湖 边界 范围 
内 的 湖面 足迹 点 , 萄 除 高 程 异常 的 湖面 足迹 点 后 计 
算得 到 湖泊 2003 一 2009 年 的 卫星 测 高 水 位 ,并 根据 
实测 水 位 数据 进行 精度 验证 ,结合 实测 数据 以 及 2 
种 测 高 数据 提取 的 水 位 信息 得 到 1970 一 2020 年 巴 
尔 喀什 湖 的 水 位 序列 ,水 位 数据 的 具体 组 成 结构 如 
图 3 所 示 。 此 外 ,采用 MOD09A1 数 据 提取 2000 一 
2020 年 湖泊 面积 数据 ,结合 相关 文献 资料 得 到 长 时 


ICESat-1 [ICESat—1 GLA14% | CryoSat— 


序 湖泊 面积 数据 ,最 后 根据 公式 ,估算 得 到 巴尔 喀 
什 湖 的 水 量变 化 。 
2.2.1 水 位 数据 提取 

(1) 卫星 测 高 数据 预 处 理 

ICESat-1 测 高 数据 的 预 处 理 主要 包括 椭 球 转 
换 和 饱和 校正 。 由 于 ICESat-1 5j CryoSat-2 M i5 
据 基 于 不 同 的 参考 椭 球 ,因此 需要 将 ICESat-1 测 高 
数据 的 参考 椭 球 (TIP 椭 球 ) 转 换 成 WGS84 椭 球 ,以 
消除 椭 球 差异 带 来 的 影响 。 此 外 ,ICESat-1 在 进 
行 高 程 测量 时 可 能 因 波 形 饱和 现象 导致 测量 足迹 
点 的 高 程 低 于 实际 高 程 值 ,因此 需要 对 其 进行 饱和 
校正 。 足 迹 点 高 程 的 计算 原理 如 公式 (1) 所 示 : 


-2 [ CryoSat-2 GDR 数 据 | , [MOD09A1 数 据 


OE CS & BOMS & EGM96 
dogs s cone] 


Baltic "Tesi" WGS84 & EGM96 "me 
& EGM2008 
m 化 水 体 
Wu 各 项 误差 改正 ERR) 指数 


i M 
提取 2003 一 2009 年 bI :提取 2000 一 2020 年 
卫星 足迹 点 高 程 bI : 9 月 湖泊 面积 


<a (NDW) 


A SLSR CES COS HP B 
明显 异常 值 ， E 
每 天 的 人 dd 


2003—2020% 
水 位 序列 


根据 湖泊 缓冲 区 提取 湖面 数据 点 


| en 
NE. XX po 准则 剔除 EL 


减 去 平均 差 值 , 变换 到 与 
实测 水 位 相同 的 水 准 面 


SERED, 
提取 平均 差 值 
Y 


1970 一 2020 年 巴尔 喀什 湖水 位 序列 


! M 
一 
| 喀什 湖面 积 序列 


巴尔 喀什 湖水 量变 化 
图 2 技术 路 线 


Fig.2 Flowchart of synthesizing multi- 


实测 数据 


source lake-level data in this study 


CryoSat-2 
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图 3 水 位 数据 组 成 结构 


Fig.3 Data composition structure diagram 
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H=h-d_deltaEllip+d_satElevCorr—N — (1) 
Xp : H SET EGM96 Kl 7k HE TRES TE SS 5h HÆF 
T/P HABER AY SFE; d deltaEllip X T/P fi Ek fil WGS84 
椭 球 之 间 的 差 值 ; d satElevCorr 为 饱和 校正 参数 ， 
可 从 数据 集 里 获得 ;NN 为 当地 大 地 水 准 面 差距 ,可 由 
MATLAB 中 的 geoidheight 函数 计算 得 到 。 

CryoSat-2 测 高 数据 的 预 处 理 主要 包括 各 项 偏 
差 改正 。CryoSat-2 卫星 雷达 脉冲 信号 在 传播 过 程 
中 发 生 的 散射 或 折射 会 影响 其 传播 速度 ,使 观测 信 
号 的 往返 时 间 产 生 延 迟 ,同时 受 各 种 自然 因素 的 影 
响 , 距 离 估 计 值 会 产生 偏差 ,因此 必须 加 以 改正 ”。 
湖泊 水 位 的 计算 公式 如 下 : 

H=H,,-R-AR-N (2) 

式 中 :为 基于 EGM96 大 地 水 准 面 的 正高 ;到 为 卫 

星 质心 到 参考 椭 球 面 的 高 度 ;R 为 卫星 到 湖面 的 距 

离 ; AR 为 各 项 误差 改正 ;N 为 当地 大 地 水 准 面 差距 ， 
可 由 MATLAB 中 的 geoidheight 函数 计算 得 到 。 

AR = Dry + Wet + Ion + Sol + Pol (3) 

式 中 :Dry 为 干 对 流 层 改 正 ; Wet 为 湿 对 流 层 改正 ; 

Ion 为 电离 层 改正 ;Sol 为 固体 潮 改正 ;Pol 为 极 测 改 正 。 

(2) 高 程 异 常 点 移 除 

数据 预 处 理 完 成 后 ,对 获得 的 所 有 水 位 值 进 
行 处 理 ,以 提取 湖泊 的 水 位 序列 ,具体 的 移 除 步 又 
Ag] eem. 

第 一 步 , 对 湖泊 边界 向 湖 心 方向 做 200m 的 组 
冲 区 ,根据 缓冲 区 筛选 卫星 足迹 点 ,以 确保 数据 点 
完全 落 入 湖 内 ,减少 可 能 与 湖岸 接触 的 高 程 点 对 单 
天 水 位 数据 的 干扰 。 

第 二 步 , 对 所 得 的 水 位 数据 点 先进 行 目 视 解 
译 , 将 与 大 多 数 水 位 值 相 差 儿 十 米 甚 至 几 百 米 偏差 
极 大 的 异常 值 剔 除 。 

第 三 步 , 用 3c 准则 剔除 单 天 水 位 数据 中 的 异 
常 值 ,然后 将 剩余 有 效 水 位 值 取 平 均 作为 日 均 水 
位 。3c 准则 的 具体 判别 方法 如 公式 (4) 所 示 : 对 于 
采集 的 数据 样本 (xx x3,…, x,) , 求 取 算数 平均 值 
x 和 剩余 误差 值 w, 即 可 求 得 均 方 根 偏差 or : 


v=x,- (4) 


P |x,-a|>30 , 则 相对 而 言 误差 较 大 , x HE 
弃 , 且 误差 大 于 3o 的 观测 数据 出 现 的 概率 为 


0.003 ; i |x,-x]<3o , 则 x 为 正常 值 ,应 保留 。 

以 ICESat-1 在 2004 年 5 月 23 日 的 过 境 点 为 例 ， 
经 30 准则 剔除 异常 点 之 后 ,高 于 343.4m 的 水 位 值 
被 剔除 ,剩余 水 位 点 在 343.17~343.39 m 之 间 波 动 
(图 4)。 


343.45 - (a) 3c 准 则 剔除 异常 值 前 


水 位 /m 


0 3 6 09 12 15 18 21 24 
时 间 /h 


343.40 - (b) 3o 准 则 剔除 异常 值 后 


时 间 /h 


图 4 3c 准则 剔除 异常 值 前 后 对 比 
Fig. 4 Comparison before and after removing outliers 


under the 3c criterion 


S E , XTA AY H IKA , Fe H a Ae El ER 
明显 离 群 值 ,再 用 30 准则 , dE — 2b Sal BR oe h E, Ba 


后 计算 得 到 月 均 水 位 与 年 均 水 位 。 

(3) 数据 转换 

如 技术 路 线 图 (图 2) 所 示 ,ICESat-1 数 据 采用 
了 TP、EGM2008 参考 系统 , CryoSat-2 数 据 采用 了 
WGS84 .EGM96 参 考 系统 ,而 水 文 站 的 实测 数据 采 
用 波罗的海 (Baltic ) 高 程 系统 ,3 种 数据 采用 了 不 同 
的 高 程 基准 ,因此 在 构建 水 位 序列 前 需要 将 其 进行 
融合 。 在 本 研究 中 ,首先 ,将 ICESat-1 数据 转换 为 
与 CryoSat-2 相 同 的 参考 系统 ,通过 上 文 的 高 程 异常 
点 移 除 方法 进行 数据 筛选 ,获得 巴尔 喀什 湖 2003 一 
2020 年 湖泊 年 均 水 位 序列 。 根 据 该 水 位 序列 与 对 
应 年 份 实测 水 位 间 的 平均 差 值 ,将 高 度 计 获 得 的 
水 位 减 去 平均 差 值 转换 为 与 实测 水 位 相同 的 水 准 
面 ,进而 得 到 巴尔 喀什 湖 1970 一 2020 年 的 水 位 
序列 。 

(4) 精度 验证 

本 文 参照 前 人 方法 ,通过 相关 系数 (7) .显著 性 
值 (P) ,平均 绝对 误差 (MAE) 以 及 均 方 根 误差 
(RMSE) 来 计算 雷达 高 度 计 提取 的 湖泊 水 位 与 对 应 
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的 实测 水 位 间 的 相关 性 。 相 关系 数 (7) 是 用 以 反 
映 变 量 之 间 相 关 关 系 密切 程度 的 统计 指标 。 


n 


2 (a-z) 


r= = (5) 
DE 

式 中 : r 为 相关 系数 ; x, 为 高 度 计 获得 的 水 位 值 ; x 
为 高 度 计 获得 的 水 位 值 的 算数 平均 值 ; y, 为 实测 水 
位 值 ; y 为 实测 水 位 值 的 算数 平均 值 ; n 为 样本 数 。 
r>0 时 表明 变量 之 间 为 正 相 关 , 反 之 为 负 相关 , r 的 
绝对 值 越 接 近 于 1, 表明 变量 间 相 关 关 系 越 强 。 

显著 性 值 (P) 是 用 以 反映 某 一 事件 发 生 的 可 能 
性 大 小 的 统计 指标 ,P<0.05 为 显著 ,P<0.01 为 非常 
显著 ,表明 样本 间 的 差异 由 抽样 误差 所 致 的 概率 小 
于 0.05 或 0.01。 

平均 绝对 误差 (MAE) 表 示 提 取 值 和 实测 值 之 
间 绝 对 误差 的 平均 值 , 均 方 根 误差 (RMSE) 是 提取 
值 与 实测 值 偏差 的 平方 与 样本 个 数 比 值 的 平方 根 ， 
2 个 指标 都 能 够 很 好 地 反映 出 提取 值 与 实测 值 之 间 
的 偏差 。 
2.2.2 面积 数据 提取 巴尔 喀什 湖 1970、1975、 
1980 , 1985 , 1990 .1995 年 面积 数据 来 源 于 相关 文 
HRE? , 2000—2020 年 逐年 遥感 影像 面积 是 基于 9 
H MODO09A1 数据 ,采用 Mcfeeters ^ FE H4 AY J — 1b 
差异 水 体 指数 (Normalized difference water in- 
dex, NDWI) 提 取 湖 泊 边 界 计算 所 得 ,NDWI 表 示 
如 下 : 


_ Green - NIR 
NDWI= Green + NIR (6) 


式 中 : Green 为 绿色 波段 , 即 MOD0O9A1 第 4 波段 ; 
NIR 为 近 红 外 波段 即 MOD09A1 第 2 波段 。 利 用 
NDWI 可 提取 水 体 信息 ,减弱 地 表土 壤 、 植 被 等 信息 
特性 的 干扰 进行 浆 值 分 割 ” ,有效 地 将 水 体 同 植被 
与 土壤 明显 区 分 开 。 
2.2.3 水 量变 化 估计 根据 内 陆 湖 的 水 量 平衡 ,内 
陆 湖 的 水 量变 化 是 由 湖区 面积 和 水 位 共同 决定 
的 。 我 们 使 用 以 下 方程 来 估算 水 量变 化 
AV=5XAHX(S,+S, + [S,xS;) (7) 
RP: AV 2 SAH EKER; AH 为 2 个 时 
期 的 湖泊 水 位 变化 ; S, A S, 为 2 个 时 期 的 湖区 
面积 。 


3 结果 与 分 析 


3.1 巴尔 喀什 湖水 位 变化 特征 

结合 ICESat-1 , CryoSat-2 2 种 卫星 测 高 数据 和 
水 文 站 实测 数据 ,得 到 了 巴尔 喀什 湖 1970—2020 4E 
湖泊 水 位 变化 情况 (图 5)。 在 研究 时 段 内 ,湖泊 水 
位 呈 先 下 降 后 上 升 趋势 。1970 一 1987 年 ,湖泊 水 位 
显著 下 降 (P<0.01) ,其 变化 率 约 为 1.23 m* (10a), 
水 位 总 体 下 降 约 2.21 m, 且 1987 年 湖泊 水 位 达到 最 
低 值 (340.64 m)。1988 一 2020 年 湖泊 水 位 总 体 呈 上 
升 趋势 (P<0.01) ,在 此 期 间 , 可 以 观察 到 多 次 较 短 的 
湖泊 水 位 振荡 ,增长 区 间 分 别 为 :1988 一 1990 年 、 
1992—1995 年 1999—2007 年 , 2009—2011 年 , 
2015 一 2020 年 ,其余 时 段 巴尔 喀什 湖水 位 呈 小 幅 下 
降 趋 势 。 其 中 1999—2007 年 湖泊 的 变化 趋势 与 
Propastin 等 .3 基于 TOPEX/Poseidon 和 jJason-1 卫星 
在 巴尔 喀什 湖 的 测 高 结果 基本 吻合 ,该 时 段 为 20 世 
纪 90 年 代 后 湖泊 水 位 的 显著 上 升 期 ,持续 上 升 时 间 
最 长 2009—2011 年 上 升 速率 最 高 ,为 1.66 m- 
(10a) ,局 部 峰值 达到 342.66 mo 
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图 5 1970 一 2020 年 巴尔 喀什 湖水 位 变化 
Fig. 5 Changes of water level in Balkhash Lake from 
1970 to 2020 


为 进一步 监测 巴尔 喀什 湖水 位 的 年 内 动态 变 
化 ,本 文 对 研究 时 段 内 各 月 的 水 位 进行 了 分 析 ( 图 
6)。 结 果 表 明 ,巴尔 喀什 湖 的 水 位 变化 具有 明显 的 
季节 性 特征 ,1 一 6 月 多 年 水 位 均值 旦 上 升 趋势 (P< 
0.01) ,上 且 6 月 达到 年 内 水 位 峰值 , 约 为 341.99 m, BG 
后 湖泊 水 位 呈 下 降 趋势 (P<0.01), 进 入 9 月 后 ,各 月 
平均 水 位 高 度 值 波动 幅度 逐渐 变 小 ,开始 趋 于 稳 
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" 平均 值 。 最 大 值 
* 最 小 值 一 中 位 线 
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342.0 上 


水 位 /mm 


341.5 


341.0 上 
340.5 上 
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图 6 1970 一 2020 年 巴尔 喀什 湖 逐 月 水 文 变化 


Fig. 6 Monthly hydrological changes of Balkhash 
Lake from 1970 to 2020 


定 , 水 位 均值 介 于 341.65~341.7 m ,说 明 巴 尔 喀 什 湖 
的 年 内 水 位 增长 主要 发 生 在 2 月 底 至 6 月 初 , 且 暖 
季 (4 一 10 月 ) 水 位 波动 比 冷 季 (11 月 一 次 年 3 月 ) 更 
剧烈 。 
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3.2 巴尔 喀什 湖面 积 变化 特征 

1970 一 2020 年 巴尔 喀什 湖面 积 的 长 期 变化 如 
图 7a 所 示 。1970 一 2020 年 ,巴尔 喀什 湖面 积 从 
19996 km? 减少 到 16641.93 km ,平均 变化 率 为 
-65.77 km *a ' (P«0.01) ,湖面 缩减 约 16.77% ,面积 
最 大 值 出 现在 1970 年 (19996 km?) ,最 小 值 出 现在 
1990 年 (16638.87 km?) 。1970 一 1990 年 ,巴尔 喀什 
湖 的 湖泊 面积 以 -159.86 km? - a~ BB IE BEC lli 28 AA 
(CP<0.01)。 进 入 21 世 纪 后 ,2000 一 2010 年 湖泊 面积 
有 所 回升 ,由 16672.89 km 增加 到 了 16896.44 km’, 
增加 约 0.013% ,面积 增加 速率 约 为 20.32 km a CP 
0.01) ,而 2010 年 之 后 ,湖泊 面积 呈 波 动 减少 的 趋 
势 ,萎缩 率 约 为 21.31 km? a (P«0.01) ,至 2020 年 ， 
巴尔 喀什 湖面 积 约 为 16641.93 km?。 
3.3 巴尔 喀什 湖水 量变 化 特征 

本 文 根 据 巴尔 喀什 湖面 积 数据 与 水 位 数据 , 计 
算得 到 巴尔 喀什 湖 多 年 水 量变 化 (图 8)。1970 一 


16950 F (b) 2000—20204£ 
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年 份 


图 7 1970 一 2020 年 巴尔 喀什 湖面 积 变 化 
Fig.7 Changes in the area of Balkhash Lake from 1970 to 2020 
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图 8 1970 一 2020 年 巴尔 喀什 湖水 量变 化 
Fig. 8 Changes of water volume in Balkhash Lake from 1970 to 2020 
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FHEWY 


2020 年 巴尔 喀什 湖水 量 损失 约 为 12.33 km’, 1970— 
1985 4E , YA BK BVA 1.84 km? .sa 的 速率 持续 下 
降 , 累积 变化 量 为 -27.55 km’, 1985—1990 FEKE 
大 幅 增 长 ,增长 量 达 10.39 km ,1995 一 2007 年 湖泊 
BK EV 1.86 km? .za 的 速率 持续 上 升 , 随 后 水 量变 
化 呈 波 动 下 降 趋 势 , 其 中 2011 一 2015 年 车 水 量 减少 
速率 最 快 (-2.65 km .a-) ,累积 减少 约 10.59 km; 
3.4 水 位 数据 精度 验证 

ICESat-1 和 CryoSat-2 卫星 的 观测 数据 与 水 位 
站 实测 数据 在 2003 一 2016 年 存在 部 分 重合 ,其 中 ， 
ICESat-1 获 取 时 间 主 要 集中 在 2 一 6 月 和 9 一 12 月 ， 
与 实测 数据 的 重合 值 有 27 个 ,CryoSat-2 为 逐 月 观 
测 , 从 2010 年 7 月 起 ,与 实测 数据 的 重合 值 个 数 为 
73。 因 此 ,为 验证 卫星 测 高 数据 的 准确 性 和 可 靠 
性 ,对 卫星 观测 数据 与 实测 水 位 数据 做 相关 性 分 
析 , 结 果 显 示 ( 图 9) ,巴尔 喀什 湖 的 ICESat-1 和 
CryoSat-2 卫星 测 高 水 位 与 水 文 站 实测 水 位 之 间 呈 
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显著 的 正 相 关 关 系 , 相 关系 数 分 别 为 0.874 和 0.805， 
显著 性 值 均 小 于 0.01, 且 ICESat-1 数 据 与 实测 数据 
相 比 平均 绝对 误差 (MAE) 及 均 方 根 误差 (RMSE ) 分 
别 为 0.206 和 0.241, CryoSat-2 数 据 与 实测 数据 相 比 
平均 绝对 误差 (MAE ) 与 均 方 根 误差 (RMSE ) 分 别 为 
0.153 和 0.199 ,结果 表明 利用 2 种 卫星 测 高 数据 来 
进行 巴尔 喀什 湖 长 时 间 序 列 的 水 位 变化 监测 具有 
一 定 的 可 行 性 。 

3.5 影响 因素 分 析 

3.5.1 和 气候 变化 巴尔 喀什 湖 流域 内 水 文 站 点 公布 
的 气象 资料 不 完整 , 仅 有 1970 一 2005 年 的 年 气温 与 
降水 量 数据 ,为 验证 CRUTS v4.05 气象 数据 集 在 巴 
尔 喀什 湖 流域 内 的 适用 性 ,本 文 对 流域 站 点 数据 与 
CRUTS v4.05 气象 数据 进行 了 相关 性 分 析 ,结果 显 
示 ( 图 10) ,气象 站 点 获得 的 气温 与 降水 量 数据 与 
CRUTS v4.05 数据 集 之 间 旦 显著 的 正 相 关 关 系 , 相 
关系 数 分 别 为 0.965 和 0.951。 表 明 CRUTS v4.05 数 
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图 9 精度 验证 


Fig. 9 Precision validation 
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图 10 CRUTS v4.05 数 据 集 精 度 验 证 
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Fig. 10 Precision validation of CRUTS v4.05 dataset 
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TRASH: 基于 多 源 遥 感 数据 的 1970 一 2020 年 巴尔 喀什 湖 动 态 监测 


据 集 可 用 于 讨论 巴尔 喀什 湖 流域 的 气候 变化 状况 。 

1970—2019 年 巴尔 喀什 湖 流域 多 年 平均 气温 
约 为 6.02 C ,流域 气温 在 研究 时 段 内 整体 呈 显 著 上 
升 趋势 ,通过 了 99% 置 信和 度 检 验 , 变 化 率 为 0.39 C- 
(103) (图 1la)。M-K 检 验 的 结果 表明 巴尔 喀什 湖 
流域 的 年 均 温 在 1988 年 发 生 突变 ,开始 出 现 增 温 趋 
势 ,2003 年 以 后 流域 的 年 均 温 上 升 趋势 显著 ,置信 
度 超过 95%( 图 1lb)。 巴 尔 喀什 湖 流 域 多 年 平均 降 
水 量 为 280.584 mm, 其 变化 率 为 0.06 mm * (10a)"', 
未 通过 90% 置 信和 度 检验 (图 11c)。M-K 检 验 结果 表 
明 巴 尔 喀什 湖 流域 的 年 降水 量 出 现 多 个 突变 点 , 降 
水 量 在 1978 一 1995 年 波动 幅度 较 大 ,而 在 1997 年 
之 后 开始 增多 ,2015 年 之 后 呈 显 著 增加 趋势 ,置信 
度 达 到 95%( 图 11d)。 

对 巴尔 喀什 湖 流 域 气温 及 降水 量变 化 与 湖泊 
动态 变化 分 时 段 进行 相关 性 分 析 ,结果 表明 水 位 、 
面积 及 水 量变 化 在 1970 一 1989 年 与 气温 降水 量 均 
无 相关 关系 ,1990 一 2019 年 , 仅 水 位 与 气温 呈正 相 
关 , 相 关 性 为 0.387( 在 95% 水 平 上 显著 )。 因 而 , 巴 
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尔 喀 什 湖 流域 气温 及 降水 量变 化 与 湖泊 水 位 .面积 
水 量 的 年 际 变化 在 长 期 趋势 上 并 无 明显 相关 性 。 

3.5.30 人 类 活动 ”1970 年 卡 普 恰 盖 水 库 建 成 后 , 伊 
犁 河 人 湖水 量 由 1950—1969 年 的 148.2x108mm 锐 减 
A| 1970—1985 4E. fff] 116.1x 10 m? ^? ,同时 随 着 流域 
灌溉 农业 的 发 展 ,湖区 及 周边 地 区 的 农业 用 水 大 幅 
增长 所 ,最 终 使 得 巴尔 喀什 湖水 位 、 面 积 以 及 水 量 
在 1970 一 1987 年 均 持 续 下 降 。1986 年 起 哈萨克 斯 
坦 有 关 部 门 终止 水 库 著 水 计划 .调节 水 库 发 电量 并 
改善 河道 淤积 状况 ,这 一 系列 措施 使 得 巴尔 喀什 湖 
生态 环境 有 所 改善 。1991 年 苏联 解体 ,哈萨克 斯 坦 
经 济 受 创 , 国 家 终止 了 对 农场 的 支持 ,实施 土地 私 
有 化 ,并 开始 征收 水 费 ,巴尔 喀什 湖 流域 的 灌溉 面 
积 因 此 直线 下 降 ”… ,许多 依靠 电力 驱动 水 泵 的 灌溉 
区 被 完全 放弃 ,农业 用 水 的 减少 进一步 促使 水 位 组 
慢 上 升 ,面积 与 水 量 随 之 呈现 出 回升 趋势 ， 
2000 一 2006 年 ,流域 灌溉 面积 介 于 大 约 2000~3000 
km 之 间 , 与 20 世 纪 90 年 代 相 比 大 幅 减 少 (图 12) 
在 此 期 间 ,巴尔 喀什 湖 的 水 位 有 了 明显 的 回升 趋 
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图 11 1970 一 2020 年 巴尔 喀什 湖 流域 气候 变化 特征 
Fig. 11 Characteristics of climate change in Balkhash Lake Basin from 1970 to 2020 
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图 12 1970 一 2006 年 巴尔 喀什 湖 流 域 灌溉 面积 


Fig. 12 Irrigated area in the Balkhash Lake Basin from 1970 to 2006 


势 。 因 而 ,巴尔 喀什 湖 流域 内 强烈 的 人 类 活动 ( 包 
括 卡 善 恰 盖 水 库 蓄 水 和 哈萨克 斯 坦 境 内 灌区 耕地 
扩张 等 ) 是 1970 一 1987 年 湖泊 动态 变化 的 主要 影响 
因素 ,同时 ,人 类 活动 也 在 湖泊 动态 变化 的 长 期 趋 
势 中 发 挥 着 重要 的 影响 作用 。 

以 往 的 研究 结果 表明 2 ,巴尔 喀什 湖 动 态 
变化 受气 温 .降水 量 .蒸发 量 冰川 融 水 以 及 人 类 活 
动 等 多 种 因素 的 综合 影响 ,定量 地 探讨 各 因素 在 不 
同时 期 湖泊 动态 变化 中 的 驱动 作用 ,系统 地 掌握 和 
理解 巴尔 喀什 湖 动态 变化 的 过 程 及 原因 ,还 需 做 更 
进一步 的 研究 。 


4 结论 

(1) 据 ICESat-1 与 CryoSat-2 的 监测 结果 ,巴尔 
喀什 湖水 位 于 1970 一 1987 年 快速 下 降 , 平 均 变化 速 
率 为 1.23 m* (10a)"',1987 年 至 今 ,巴尔 喀什 湖水 位 
总 体 呈 现 上 升 趋势 ,并 伴随 有 周期 性 波动 。 同 时 水 
位 具有 较为 规律 的 季节 性 变化 ,春季 多 年 平均 水 位 
最 高 ,秋季 最 低 ,水 位 增长 主要 发 生 在 湖泊 结 冰 期 
(11 月 底 到 次 年 4 月 初 ), 暧 季 (4 一 10 月 ) 水 位 波动 
比 冷 季 (11 月 一 次 年 3 月 ) 更 剧烈 ;研究 时 段 内 ,巴尔 
喀什 湖面 积 呈现 出 先 明 显 萎 缩 后 略 有 恢复 的 变化 
过 程 ,51 a 湖面 缩减 约 16.77%; 在 此 期 间 , 巴 尔 喀 什 
湖水 量 损 失 约 为 12.33 km’, 

(2) 巴尔 喀什 湖 动态 变化 受 


气候 变化 与 人 类 活 


动 的 共同 影响 ,1970 一 1987 年 水 位 .面积 及 水 量变 
化 主要 由 卡 普 恰 盖 水 库 蔷 水 和 哈萨克 斯 坦 境 内 灌 
区 耕地 扩张 等 人 类 活动 引起 ,但 研究 时 段 内 ,巴尔 
喀什 湖 的 动态 变化 与 气候 变化 在 总 体 趋势 上 无 明 
显 相 关 性 ,定量 描述 二 者 之 间 的 联系 ,还 需 进一步 
探讨 。 

(3) 测 高 卫星 为 湖泊 水 位 监测 提供 了 有 效 的 技 
术 手 段 ,与 常规 水 位 测量 方法 相 比 , 测 高 卫星 不 受 
自然 环境 及 人 力 因 素 的 限制 ,可 有 履 盖 更 多 湖泊 , 整 
合 不 同 卫星 的 测 高 数据 ,还 可 获得 较 长 时 序 的 湖泊 
水 位 信息 ,这 对 研究 湖泊 水 位 的 动态 变化 及 其 驱动 
因素 具有 重要 意义 。 
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Abstract: Due to climate change and anthropogenic activities, water resource issues are challenging and have 
attracted significant attention worldwide. Lake level variation is an essential indicator of regional and global 
climate and environment. For the past few decades, the development of remote sensing has enabled long-term 
dynamic monitoring of lake levels and areas. In this work, the water level change of Balkhash Lake during the 
period of 1970—2020 was derived from ICESat-1, CryoSat-2, and the hydro-meteorological stations. Moreover, 
combined with MODO9A1 data and related literature, the area and water volume change of Balkhash Lake in a 
long time series were analyzed. The results indicated that during the study period, the water level of Lake 
Balkhash first showed a downward trend, and reached its lowest point (340.64 m) in 1987, and then began to 
fluctuate upward. Generally, the water level growth during the year mainly occurred from the end of February to 
the beginning of June. The water level fluctuates more dramatically in the warm season (from April to October) 
than in the cold season (from November to March of the following year). The multi-year average water level is 
the highest in spring and lowest in autumn. From 1970 to 2020, the area of the Balkhash Lake decreased from 
19996 km’ to 16641.93 km’, and the lake area decreased by about 16.77%. Meanwhile, the water volume of the 
Balkhash Lake undergone a total loss of approximately 12.33 km’. The analysis revealed that the dynamic 
changes of the lake from 1970 to 1987 were mainly caused by human activities, such as the impoundment of the 
Kapchagay Reservoir and the expansion of arable land in the irrigation area in Kazakhstan. However, the overall 
trend of temperature and precipitation has no obvious correlation with the dynamic changes of the lake. Owing to 
the complexity and uncertainty of environmental changes, quantitative discussion of the impact of various factors 
on lake dynamics requires further research. Multi-mission altimetry satellite provides a powerful tool for the long- 
term and large-scale monitoring of lake level variations, which is of great significance to the study of lake level 
change and its response to climate and environment. 
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